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ABSTRACT 

The capability of the reconstruction scheme developed in Part I is demonstrated here through three practical 
applications. First, the nonlinear regression model is used to reproduce the upper-layer three-dimensional cir¬ 
culation of the eastern Black Sea from model data distorted by white and red noises. Second, the quasigeostrophic 
approximation is used to reconstruct the shallow water circulation pattern in an open domain with various 
sampling strategies. Third, the large-scale circulation in the Southern Ocean is reproduced from the First Global 
Atmospheric Research Program (GARP) Global Experiment (FGGE) drifter data with noncontrollable noise 
statistics. All three cases confirm that the theoretical approach is robust to various noise-to-signal ratios, number 
of observations, and station disposition. Using the simplified open boundary condition for analyzing long-term 
observational data is recommended because the nonlinear regression procedure requires considerable computer 
resources. 


1. Introduction 

In Part I (Chu et al. 2003), a theoretical approach was 
developed for reconstructing a three-dimensional ve¬ 
locity field from sparse and noisy data. The scheme can 
be used for reconstructing velocity field from data in 
Eulerian and Lagrangian frames. To demonstrate the 
usefulness of this approach and to test the sensitivity of 
the approach to possible variations of noise-to-signal 
ratio, station disposition, and observational sample size, 
the schemes are applied to reconstruct the circulation 
in the eastern Black Sea, a rectangular open domain, 
and the Southern Ocean. The first two cases (identical 
twin experiments) represent the reconstruction of ve¬ 
locity from Eulerian “observations,” and the third one 
(using observational data) represents the reconstruction 
of velocity from Lagrangian data. 

The first experiment is designed to verify the model 
sensitivity to the selection of optimal model parameters 
( K opt , M opt , S opt ) through reconstructing the circulation 
in the eastern Black Sea. It is of considerable interest 
because the pipeline “Blue Flow” between Russia and 
Turkey will be built in the near future. Assessing the 
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ecological safety of the project has become urgent, and 
a realistic three-dimensional flow field is needed. Re¬ 
construction of the three-dimensional baroclinic circu¬ 
lation in the eastern part of the Black Sea was performed 
with realistic topography and controlled noise in the 
observations. The spatial structure of the circulation was 
obtained using both scalar potentials M? and O. 

The second experiment is designed to verify the mod¬ 
el sensitivity to the sampling strategy through recon¬ 
structing circulation in a rectangular open domain with 
a flat bottom. Such an experiment has a practical sig¬ 
nificance. Presently, high-frequency radars are widely 
used in measuring circulation in the coastal ocean. Usu¬ 
ally, the radar data are collected in a domain of limited 
size with high-resolution grids, which is very attractive 
if the circulation can be reconstructed in a domain larger 
than the radar polygon (Shulman et al. 2000). The ques¬ 
tion arises, how large should the domain be for a realistic 
flow reconstruction? In the experiment, we controlled 
the noise level in the data and used a simplified con¬ 
dition at the open boundary. Preliminary estimations 
demonstrated that, with given wind forcing, subgrid pa¬ 
rameterization, and basin geometry, the shallow water 
circulation can be described using the quasigeostrophic 
streamfunction, which coincides with the toroidal po¬ 
tential TL 

The third experiment is designed to verify the scheme 
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Fig. 1. The exact velocity field at 20 m deep in the eastern Black 
Sea. Note that T[ and T' 2 are the open boundary and extended open 
boundary, respectively. 


using observational buoy data with unknown noise sta¬ 
tistics and to examine the validity of boundary extension 
method through reconstructing the large-scale circula¬ 
tion in the Southern Ocean. More than 200 drifting 
buoys were deployed in the Southern Ocean during the 
First Global Atmospheric Research Program (GARP) 
Global Experiment (FGGE) from December 1978 to 
December 1979. Such a dataset provides useful infor¬ 
mation on the seasonal variability of circulation from 
70° to 20°S (Hoffman 1985; Patterson 1985; Daniault 
and Menard 1985; Eremeev and Ivanov 1987). Here, 
we use the FGGE buoy data to reconstruct mean surface 
circulation in the Southern Ocean from 70° to 48°S on 
the 2° X 2° grids in May 1979. The oceanographic in¬ 
terpretation of the reconstructed seasonal variability of 
the circulation will be presented in another paper. 

2. Case 1: Eastern Black Sea circulation 



Fig. 2. Four station dispositions: (a) 22 stations with data at the 
open boundary, (b) 19 stations without data at the open boundary, 
(c) 14 stations without data at the open boundary, and (d) 9 stations 
without data at the open boundary. 


water mass and with assimilation of climatological 
monthly mean three-dimensional temperature and salin¬ 
ity data (Eremev et al. 1992; Staneva et al. 1995); see 
Knysh et al. (2001) for detailed information about the 
model and numerical integration. 

The final solutions show a quasiperiodic characteristic 
in the circulation pattern. We use the three-dimensional 
circulation calculated at 20 m deep on 15 December as 
the exact “observational” velocity field, which shows 
a cyclonic gyre centered near 42°N, 39°E in the eastern 
Black Sea with a jet stream (rim current) and several 
small anticyclonic eddies (around 100 km in diameter) 
around it (Fig. 1). This simulated “observational da¬ 
taset,” treated as “perfect” (no error), is used as a data 
source. 


The first experiment is designed to determine the 
model sensitivity to the selection of optimal model pa¬ 
rameters (K opt , M opt , 5 opt ) and to verify the boundary 
extension method through reconstructing the circulation 
in the eastern Black Sea. 


a. Velocity “data” 

A z-coordinate primitive model is used for the whole 
Black Sea with 14.6-km horizontal resolution and 21 
vertical levels (Knysh et al. 2001) to obtain the three- 
dimensional circulation. The horizontal and vertical vis¬ 
cosities are taken as 10 7 cm 2 s _1 and 5 cm 2 s -1 for the 
entire water column. The horizontal diffusion coefficient 
is set to be 5 X 10 5 cm 2 s -1 , and the vertical diffusion 
coefficient varies with depth in the following way: 2.5 
cm 2 s _1 above 5 m deep, 0.1 cm 2 s _1 between 5 and 
15 m deep, and 0.03 cm 2 s _1 below 15 m deep. The 
model is enclosed by rigid lateral and bottom boundaries 
with no-slip and no heat and salt exchange conditions. 
The model is integrated for 25 yr from the rest with 
monthly mean surface fluxes of momentum, heat, and 


b. Controlled noise 


Four different sampling strategies A, B, C, and D 
(Fig. 2) are used to obtain data from a velocity data 
source (Fig. 1). Notice that there are no “observations” 
at the open boundary T[ for B-D data sampling. The 
station dispositions are designed to investigate the effect 
of sampling distribution and density on the reconstruc¬ 
tion quality. 

An error (noise) is added to the sampled observational 
data. The noise is represented by a two-dimensional, 
isotropic, incompressible, and homogeneous Gaussian- 
type field with the correlation function 


G(r - r') = / 2 exp 


(r - r') 2 
T 2 

^cor 


( 2 . 1 ) 


Here, r = xi + yj is the displacement vector; L cor and 
I are the correlation radius and intensity of noise, re¬ 
spectively; L cor —» 0 for the white noise and otherwise 
for the red noise. The stochastic field with the corre- 
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Fig. 3. Typical observational data collected from (a) disposition A 
with a white noise ( rj 1 = 0.3), (b) disposition A with a color noise 
( r] 1 = 0.3, L cor = 45 km), (c) disposition A with a strong white noise 
(r] l = 1.2), (d) disposition B with a white noise ( rj l = 0.3), and (e) 
disposition C with a white noise ( r] l = 0.3). 


lation function (2.1) was used in Sabelfeld’s (1991) tech¬ 
nique. The observations are distorted by various noises. 

Three typical sets of sampled data (A, B, C) with 
white and red noises are considered (Fig. 3). Notice that 
if high-intensity noise is introduced into the perfect ob¬ 
servations, the reconstruction process becomes weakly 
sensitive to whether there are observations along the 
open boundary or not. 


c. Noise-to-signal ratio 


Two types of error affect the noise-to-signal ratio : 
stochastic noise in the distorted observations, and mode 
truncation [see (3.3) in Part I]. The second type of er¬ 
ror—that is, cutoff of high-order modes—is regarded as 
the deterministic perturbation. 

The noise-to-signal ratio is represented by (Eykhoff 
1973): 


Vi = 


I+VE^ t 

Ve 


( 2 . 2 ) 


where E and E cut are the total and truncated horizontal 
circulation energies, respectively. In the numerical ex¬ 


periments r] l varies from 0.1 to 1.2. When the toroidal 
(poloidal) potential is truncated at around 20 (3), the 
contribution of the deterministic perturbation to may 
be neglected. 

d. Toroidal and poloidal components 

In the reconstruction, the open boundary extension 
method was used (Fig. 1). With boundary conditions 
given along T ' 2 , 200 basis functions are calculated for 
both scalar potentials using MATLAB (1997). The nu¬ 
merical procedure for constructing the rotation matrix 
R is discussed by Ivanov et al. (2001). The integral 
equation (3.5) of Part I is reduced to linear algebraic 
equations and solved using a package of computational 
programs from Tikhonov et al (1990). Here, this ap¬ 
proach is compared with the smoothed splines only. 

For scalar ocean fields, the regularization technique 
in a fitting model is compared with the optimal inter¬ 
polation (OI), the spline interpolation based on cross 
validation (CV), generalized cross validation (GCV), 
and Tikhonov’s approach (Ivanov et al. 2001). Shultz 
(1999) compared the reconstruction technique based on 
the mode decomposition (Eremeev et al. 1992) with OI 
for reconstructing surface currents on the Louisiana- 
Texas shelf. 


e. Optimal truncation 

A question that needs to be answered before recon¬ 
struction is how to select the optimal model parameters: 
K opt and M opt , the number of toroidal and poloidal 
modes; and S opt , the order of continuous polynomial 
used to parameterize the scalar function k(t). Vapnik- 
Chervonkis cost function (VCCF) is used to determine 
(K opt , M opt , S opt ) as explained in section 5a of Part I. For 
disposition A and B, the ranges of truncation parameters 
are found: 

K opi = 15 - 23, M opt =1-4, and S opt = 5, 

which does not vary with the data noise level noticeably. 
For simplicity, the selection of K opt and M opt is dem¬ 
onstrated using the empirical (7 emp ), VCCF (J NM ), and 
statistical (J) cost functions for a given S opt = 5, with 
the observational data distorted by 30% white noise 
(Fig. 3a). 

1) Empirical cost function 

Figure 4a shows the dependence of the empirical 
functional / emp on the mode truncation parameters K and 
M. A maximum value of J emp is found at (^T opt , M opt ) = 
(19, 4). Using 

(K opt , M opt , S opt ) = (19, 4, 5), 

the reconstructed field has an unrealistic jet near Cau¬ 
casus, which indicates the overdetermination of the po¬ 
loidal component in the general circulation pattern. 
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2) VCCF 

The dependence of VCCF on the mode truncation 
parameters K and M is demonstrated in Fig. 4b. Two 
maximum values of VCCF are found: 

(K opt , M opt \ = (17, 1), and (K opt , M opt ) 2 = (22, 1). 

It is also noticed that the mode truncations varying from 
CK opt , M opt \ to (^ opt , M opt ) 2 does not change the VCCF 
noticeably. 

3) Statistical cost function (J) 

Figure 4c shows the dependence of the statistical cost 
function (J) on the mode truncation parameters K and 
M. Three maximum values of (J) are found: 

(K opt , M opt \ = (17, 2), (^ opt , M opt ) 2 = (22,1), and 

(K opt , M opt ) 3 = (23, 0). 

The statistical cost function (J) reaches the global ex¬ 
tremum for (^ opt , M opt ) 3 . Notice that (J) is constructed 
from 1000 stochastic realizations of the distorted ob¬ 


servations and, additionally, averaged over all grid 
points used for numerical circulation modeling. 

f Verification of reconstruction scheme for 
disposition A 

We use the truncation 

K opt = 21, M opt = 1, and S opt = 5 (2.3) 

to reconstruct the observational data with various sto¬ 
chastic noises. 

1) White noise (r] l =0.3) 

A white noise with = 0.3 is added to the data 
collected from the disposition A (Fig. 2a). A typical 
distorted observational dataset for the reconstruction is 
illustrated in Fig. 3a. 

(i) Zero-order approximation h [0] 

Using the zero-order approximation h [0] , the recon¬ 
structed circulation (Fig. 5a) has a similar pattern as the 
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Fig. 5. Reconstructed velocity field from observed field indicated in Fig. 3a using the zero-order approximation: (a) velocity vectors, (b) 
residue velocity vectors, (c) normal velocity, and (d) tangential velocity at the open boundary. Here, the cross (+) and solid curve represent 
the reconstructed and exact velocity components at the open boundary. 


exact field (Fig. 1): an areawide cyclonic gyre centered 
at 41.2°N, 39.15°E with a jet (rim current) around it. 
The center of the cyclone is shifted 20 km toward the 
southeast. The residue (exact minus estimated) velocity 
vector field (Fig. 5b) shows that the scheme overesti¬ 
mates the eastward jet near the southern boundary of 
the Black Sea and underestimates the westward jet in 
the northwestern part of the region. Besides, an anti¬ 
cyclone near the eastern corner of the Black Sea (cen¬ 
tered at 42.7°N, 41.0°E, and 100 km in diameter) oc¬ 
curring in the reconstruct field (Fig. 5a) does not exist 
in the exact field (Fig. 1). 

The capability of the reconstruction scheme is also 
recognized by the agreement between the reconstructed 
(cross) and exact (solid) normal (Fig. 5c) and tangential 
(Fig. 5d) velocity components at the open boundary. 
The reconstruction skill is represented by the root-mean- 
square (rms) errors [see (5.10) and (5.11) of Part I]: 

Xu = 0.13, Xv = 0.18, Xn = 0.08, and 

Xr = 0.25. (2.4) 

Hence, the reconstruction quality is much higher for the 
normal component than for the tangential component. 


(ii) Nonlinear regression 

Using the nonlinear regression (see section 4b in Part 
I) with optimal truncations (2.3), the reconstructed cir¬ 
culation (Fig. 6a) has a pattern closer to the exact ob¬ 
servational field (Fig. 1) than the zero-order approxi¬ 
mation: a basin wide cyclonic eddy centered at 41°N, 
39°E with a jet (rim current) and several small eddies 
around it. The center of this cyclone has the same lo¬ 
cation as the exact observational field. The residue ve¬ 
locity vector field (Fig. 6b) shows that the scheme un¬ 
derestimates the westward jet in the northwestern part 
of the region. Similar to the zero-order approximation, 
an anticyclone near the eastern corner of the Black Sea 
(centered at 42.5°N, 41.2°E, and 100 km in diameter) 
occurs in the reconstruct field (Fig. 6a), which does not 
exist in the exact observational field (Fig. 2). However, 
the strength of this false anticyclone is weaker than that 
in the zero-order approximation. 

The reconstructed (cross) normal (Fig. 6c) and tan¬ 
gential (Fig. 6d) components agree well with the exact 
(solid) field at the open boundary. The reconstruction 
skill is represented by 
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Xu = 0.10, Xv = 0.13, Xn = 0.11, and 

X, = 0.14, (2.5) 

which indicates that the rms errors of all velocity com¬ 
ponents are less than or equal to 0.14. Comparison be¬ 
tween (2.4) and (2.5) shows the improvement in recon¬ 
struction (except the normal velocity at the open bound¬ 
ary) using the nonlinear regression. 

From theoretical point of view, the nonlinear regres¬ 
sion should provide more accurate results than the zero- 
order approximation. However, comparison between 
(2.4) and (2.5) shows more accuracy on the normal ve¬ 
locity at the open boundary using the zero-order ap¬ 
proximation rather than the nonlinear regression. This 
may be caused by strong dependence of the reconstruc¬ 
tion accuracy on various factors such as station depo¬ 
sition pattern, level of noise in observations, and quality 
of numerical algorithm; further investigation is needed. 

2) Color noise 

A color noise with = 0.3, L cor = 45 km is added 
to the data acquired from the disposition A (Fig. 2a). 
A typical distorted observational dataset for the recon¬ 
struction is illustrated in Fig. 3b. 


(i) Zero-order approximation 

Using the zero-order approximation h [0] , the recon¬ 
structed circulation (Fig. 7a) has a pattern similar to the 
exact circulation field (Fig. 1): a basinwide cyclonic 
gyre centered near 42°N, 39°E with a jet stream and 
several small eddies (about 100 km in diameter) around 
it, and this cyclone seems to be flattened. The east-west 
semiaxis is much longer than the north-south semiaxis 
(Fig. 7a). The residue velocity vector field (Fig. 7b) 
shows that the scheme overestimates both the eastward 
jet near the southern boundary and the westward jet in 
the northwestern part of the region. Note that an anti¬ 
cyclone near the upper half of the eastern corner of the 
Black Sea (centered at 42.75°N, 41°E, and 100 km in 
diameter) occurring in the reconstruct field (Fig. 7a) 
does not exist in the exact observational field (Fig. 1). 
At the open boundary, the reconstruction quality for the 
normal component (Fig. 7c) is much better than the 
tangential component (Fig. 7d). The reconstruction skill 
is represented by 

Xu = 0.17, Xv = 0-16, Xn = 0.05, and 

X, = 0.48, (2.6) 

which indicates that the reconstruction quality for the 
data with red noise is comparable to that for the data 
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Fig. 7. Same as in Fig. 5, except from observed field indicated in Fig. 3b. 


with white noise at the same noise-to-signal level (0.3) 
except for the tangential velocity at the open boundary. 
The tangential (normal) velocity error is much higher 
(lower) for the red noise than for the white noise. 

(ii) Nonlinear regression 

Using the nonlinear regression with optimal trunca¬ 
tions (2.3), the reconstructed circulation (Fig. 8a) shows 
a similar pattern to the exact observational field (Fig. 
1). The residue velocity vector field (Fig. 8b) shows that 
the scheme overestimates the eastward jet near the 
southern boundary and underestimates the westward jet 
in the northwestern part of the region. Similar to the 
zero-order approximation, an anticyclone near the east¬ 
ern corner of the Black Sea (centered at 42.6°N, 41.2°E, 
and 100 km in diameter) occurs in the reconstruct field 
(Fig. 8a), which does not exist in the exact observational 
field (Fig. 1). At the open boundary, the reconstruction 
quality for the normal component (Fig. 8c) is compa¬ 
rable to that for the tangential component (Fig. 8d). The 
rms errors are computed by 

Xu = 0.1, Xv = 0.29, Xn = 0.15, and 

Xt = 0.12. (2.7) 

As with the white noise, the nonlinear regression has 
less (more) accuracy than the zero-order approximation 


in reconstructing the normal (tangential) component at 
the open boundary. 

3) Strong white noise 

To test the reconstruction capability in processing 
highly noisy data, a strong white noise (rj 1 = 1.2) is 
added to the data acquired from the disposition A (Fig. 
2a). The distorted observational dataset for the recon¬ 
struction is illustrated in Fig. 3c. The nonlinear regres¬ 
sion with the optimal truncation (2.3) is used to recon¬ 
struct the velocity field. 

The reconstructed circulation (Fig. 9a) still shows the 
existence of a basinwide cyclonic circulation. The major 
cyclone centered near 42°N, 39°E is also reconstructed. 
The residual velocity vector field (Fig. 9b) shows that 
the scheme overestimates the eastward jet near the 
southern boundary. The exact (solid) and reconstructed 
(cross) normal (Fig. 9c) and tangential (Fig. 9d) veloc¬ 
ities at the open boundary show the capability of the 
scheme in reconstructing the velocity field at the open 
boundary. The reconstruction skill is represented by 

Xu = 0.83, Xv = 1-22, Xn = 0.91, and 
= 0.83. 


( 2 . 8 ) 
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Fig. 8. Same as in Fig. 5, except from observed field indicated in Fig. 3b and using the nonlinear regression. 


g. Sensitivity of reconstruction to sampling density 

To test the dependence of the reconstruction capa¬ 
bility on the sampling density, three observational da¬ 
tasets are acquired from dispositions B, C, and D, dis¬ 
torted by a white noise ( r] l = 0.3) (Fig. 3). The nonlinear 
regression is used to reconstruct the velocity fields. The 
rms errors (Table 1) show that the reconstruction quality 
deteriorates with the reduction of sampling density from 
high sampling density (disposition B, 19 observations) 
to medium sampling density (disposition C, 14 obser¬ 
vations), and from medium sampling density to low 
sampling density (disposition D, 9 observations). For 
the low sampling density (disposition D), only six to¬ 
roidal modes and one poloidal mode can be estimated. 
Any combination of these modes results in a poor re¬ 
construction skill. 

h. Reconstruction of vertical velocity 

1) Vertical velocity “data” 

If surface gravity waves are ignored, the surface (z 
= 0) boundary condition may be written as 

w = 0, at z = 0. (2.9) 

The “exact” vertical velocity data (w ex ) at 20 m deep 
(Fig. 10a) are obtained for verification of the recon¬ 


struction scheme from vertical integration of the con¬ 
tinuity equation of the Black Sea circulation model from 
the surface to 20 m. 

2) Two-step reconstruction 

The reconstruction process includes two steps. First, 
the observational horizontal velocity field is acquired 
from disposition A and distorted by a white noise ( r] l 
= 0.3) at z = 0, 5, 10, and 20 m, and the horizontal 
velocity is reconstructed at each level. Then the vertical 
gradient of poloidal potential (0 to 20 m) is determined. 
Notice that the vertical gradient of O does not equal 
zero at the ocean surface. Second, the integral equation 
(3.5) of Part I is solved to obtain the reconstructed ver¬ 
tical velocity w rec at 20 m (Fig. 10b). Here, strong cor¬ 
relation is found between the vertical velocity (upwell- 
ing and downwelling zones; Fig. 10b) and the eddy 
structures in horizontal velocity vector field (Figs. 5a 
and 6a). The reconstruction error estimated by 


Large values of w res (Fig. 10c) are found in some areas, 
because the reconstruction skill of the vertical velocity 
strongly depends on that of the horizontal velocity. At 
least five or six poloidal modes are needed to best re¬ 
solve the horizontal velocity structure during recon- 
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Fig. 9. Same as in Fig. 5, except from observed field indicated in Fig. 3c and using the nonlinear regression. 


struction if the observational data contains no noise. For 
the selected station disposition and noise level, a max¬ 
imum of two or three poloidal modes can be recon¬ 
structed. This limitation causes large difference between 
exact and reconstructed vertical velocity in some areas. 

i. Comparison with the smoothed splines 

The reconstructed field is compared to that repro¬ 
duced using the smoothed splines whose smoothing pa¬ 
rameter is obtained using Tikhonov’s regularization pro¬ 
cedure. Both techniques weakly depend on the noise 
level in observations if rj l < 0.3. For example, the re¬ 
construction skill of the smoothed splines over the ob¬ 
servational data acquired from disposition A and dis¬ 
torted by a white noise (= 0.3), is given by 

Xu = 0.15, Xv = 0-23, Xn = 0-18, and 

Xr = 0.43, (2.10) 


Table 1. Reconstruction skill for various disposition patterns. 


Sampling pattern 


X v 

X n 


Disposition B 

0.30 

0.23 

0.18 

0.43 

Disposition C 

0.28 

0.36 

0.49 

0.48 

Disposition D 

0.32 

0.40 

0.55 

0.62 


which is comparable to (2.4) and (2.5). For the obser¬ 
vational data collected from disposition A and distorted 
by a strong white noise (= 1.2), the rms errors are 
extremely high using the smooth splines: 

Xu = 1-18, Xv = 1-73, Xn = 0.65, and 

= 1 . 1 . ( 2 - 11 ) 

Comparison between (2.11) and (2.8) confirms the ca¬ 
pability of the reconstruction scheme using the two po¬ 
tential scalar representation. 

3. Case 2: Shallow water flow in a semienclosed 
basin 

a. Geometry of the basin 

A shallow water basin with a flat bottom is centered 
at 35°N and bounded by three rigid boundaries. This 
basin extends 1000 km in the north-south and 1050 km 
in the east-west directions. The northern, southern, and 
western boundaries are rigid, and the eastern boundary 
is open. The Cartesian coordinate system is chosen with 
the origin at the southwest corner. The v axis points 
toward the east, and the y axis toward the north (Fig. 
11a). 
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Fig. 10. Vertical velocity (in 10 -4 cm s -1 ) at 20 m deep: (a) exact 
data, (b) reconstructed field from depostion A sampling with a white 
noise ( rj 1 = 0.3), and (c) error (w res ). The positive values indicate 
up welling. 


b. Exact velocity data 

The exact velocity data for the semienclosed basin is 
model generated with the Princeton Ocean Model 
(POM) developed by Blumberg and Mellor (1987). 


(a) 



Fig. 11. Two areas for the POM integration: (a) domain with three 
rigid and one open boundaries, and (b) a domain with four rigid 
boundaries. 


POM is a primitive equation model with a free surface 
and a level 2 1/2 turbulence closure scheme. We use 
the 2D version of POM to get the exact data (Chu et 
al. 1997) for verifying the reconstruction scheme. The 
eastern boundary is connected to a mirror image of do¬ 
main A (about x = 1050 km), forming a closed rect¬ 
angular domain (Fig. lib). The POM model was inte¬ 
grated for the large domain with four rigid boundaries 
from the following initial conditions: 

(u, v, w) = 0 (in m s -1 ), 


T = 


283 


1 + exp 



(in K), 


and 


S = 35 (in ppt), 


and with no surface heat or salinity fluxes but zonal 
surface pseudo-wind stress varying with latitude: 


T 77V 

— = -10- 3 cos—, L = 1000 km. (3.1) 
Po C 

The model time step was 2.5 min. The horizontal 
resolution is 50 km. Bottom stress is parameterized by 
the quadratic drag relation and the horizontal kinematic 
viscosity is set to 0. When POM is integrated from the 
rest (u = v = 0) and the flat sea surface (£ = 0), two 
circulation regimes are identified in the basin: (i) a linear 
regime (up to day 20) with velocities under 10 to 15 
cm s _1 and surface sea elevation reaching 5-10 cm; (ii) 
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(a) 


(b) 


Circulation in channel at 5day Circulation in channel at 90day 




Fig. 12. Exact velocity dataset for the semienclosed basin: (a) day 5 and (b) day 90. 


a nonlinear regime (from day 20 to day 50) with oc¬ 
currence of large velocities (100 to 120 cm s _1 ) and 
surface elevation (100 cm). For the nonlinear regime, 
we have 


= max l jrj ~ °' 25> 


- max| -Jjr \ - 0.2, 


where U ± is the characteristic velocity in the basin; L ± 
and T ± are the spatial and temporal scales of the flow, 
respectively; / is the Coriolis parameter. The POM- 
generated 2D velocity fields on day 5 (in the linear 
regime; Fig. 12a) and day 90 (in the nonlinear regime; 
Fig. 12b) are treated as the exact velocity data. 


__i_ (<w ao 


H + £;\dy 

1 


dx / 


H+ g 


d ( t> \ 

-1-, and 

dx dy I 


— = -A<2>, 
dt 


(3.2) 


and the streamfunction 1 P satisfies the Poisson equation 

AV = %-^v - (H + (3-3) 

dy dx 


c. Sampling strategies 


where f is vertical vorticity. 
Scale analysis leads to 


The following six sampling strategies (Fig. 13) are 
used to investigate the effect of data distribution on the 
quality of reconstruction: (a) irregular basinwide dis¬ 
position with high density (86 points), (b) irregular bas¬ 
inwide disposition with lower density (42 points), (c) 
regular basin wide disposition (192 points), (d) irregular 
partial domain (noncoastal) disposition (69 points), (e) 
regular coastal disposition (114 points), and (f) irregular 
coastal disposition (69 points). Notice that there are no 
observations at the open boundary. Here, we are par¬ 
ticularly interested in the reconstruction quality for par¬ 
tial sampling cases (i.e., dispositions D, E, and F), where 
the observational domain (, S ) is smaller than the recon¬ 
struction domain ( S '), that is S C S'. 


d. Simplification of two-scalar representations 

For a 2D shallow water flow, the toroidal and poloidal 
potentials coincide with the streamfunction and velocity 
potential (Chu 1999; Chu et al. 2003). The velocity ( u , 
v) and surface elevation (£) are represented by 


if! 

H 


o(io- 4 ), 


I udg/dy — vd^/dx I 

I (H + && 


O(10- 4 ), 


\d<&/dy\ 
dV/dx I 


O(10“ 6 ), 


\d®/dx\ 

\dWdy\ 


O(10“ 6 ). 


Thus, Eq. (3.2) is simplified into 


1 d'P _ 1 dT' 

H dy 9 V ~ H dx' 


(3.4) 


Thus, the reconstruction is reduced to the analysis of 
one scalar potential (i.e., the streamfunction TO- Notice 
that the neglect of the second scalar potential (<f>) is 
accurate only for the selected rectangular basin char¬ 
acteristics. The exact data illustrated in Figs. 12a and 
12b show that at the open boundary the normal com¬ 
ponent is much larger than the tangential component | L 
I I dr » I V T I Jr at open boundary. We may ne¬ 
glect V T and use k = 0 at the open boundary [see (2.10) 
of Part I] for the reconstruction. 
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DISPOSITION A DISPOSITION B DISPOSITION C 



DISPOSITION D DISPOSITION E DISPOSITION F 



x (km) x (km) x (km) 


Fig. 13. Several sampling strategies: (a) irregular domainwide disposition with high density (86 points), 
(b) irregular domainwide disposition with lower density (42 points), (c) regular partial domainwide disposition 
(192 points), (d) irregular partial domain (noncoastal) disposition (69 points), (e) regular coastal disposition 
(114 points), and (f) irregular coastal disposition (69 points). 


e. Noise reduction using the Gram-Schmidt process 

The basis functions {^ k I k = 0, 1, . . . , K} are cal¬ 
culated using MATLAB (1997), and estimate the spec¬ 
tral coefficients a k using the linear regression model. 
The optimal mode truncation is obtained using VCCF 
[see (3.3) in Part I]. The truncation error 

'Ppernb = 2 a t % 

K 

is defined as the deterministic perturbation. The noise- 
to-signal ratio associated with the deterministic pertur¬ 
bation can reach 0.1. Such a level of noise may con¬ 
siderably deteriorate the reconstruction skill, especially 
for a partially sampled region such as dispositions C, 
D, and E (Fig. 13). 

The regularization procedure [see (4.6)-(4.9) of Part 
I] effectively filters noise if vectors Y noise and Y exact are 
perpendicular: 

Y Y =0 (3 5) 

Therefore, the reconstruction procedure is transformed 
to satisfy (3.5), which can be done through selecting 
appropriate basis functions {T^}. 

Since the number of observations (Fig. 13) is suffi¬ 
ciently large, two sets of basis functions can be con¬ 
structed from the gradients of 


© 


k 


Wk 

dy 


and 


ay, 

dx 


which can be orthonormalized using the Gram-Schmidt 
process within S. Thus, we represent each of the velocity 
components as 


K K 

U K ~ 2 a k®k an d V K — ~ 2 bk^k- 

k =0 k=Q 


(3.6) 


Next, the two sets {©^} and {H^} are reorthonor- 
malized within the domain S' into 

(H&' k ©;j = 8 km , (HZ' a'J = 8 km , 

where &' m and a' m are renormalized basis functions, 


f 1, X E S 

H = \ 

[ 0 , x £ S 

is the weight, and 8 nm is the Kroneker delta. The velocity 
components can be expanded into 

K K 

u N — 2 A m & m v N — — 2 (3.7) 

m= 0 m= 0 


and new spectral coefficients A m and B m are computed by 

K K 

K = S t^Sk a k and B m = ^j (3.8) 

k =0 k =0 











































504 


JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY 


Volume 20 


where /jl^I and /zj$ are the coefficients generated using 
the Gram-Schmidt process. Using new basis functions 
{©^}, {E'J, the coefficients A m and B m are estimated 
from the observations inside S by the simplest regression 
model. The accuracy of the reconstruction should be 
improved since the influence of deterministic pertur¬ 
bation on the coefficients has been reduced. Then, the 
spectral coefficients a k and b k are calculated from (3.8), 
which gives an approximate solution for the domain S '. 


f Reconstructed velocity field 

Since the semienclosed basin (Fig. 11) is rectangular 
and k = 0 is used for the open boundary condition [see 
(2.10) of Part I], the basis functions T ' k may be replaced 
by the trigonometric polynomials. Such a replacement 
may simplify the theoretical analysis; however, the con¬ 
vergence of the spectral expansion of the streamfunction 
is much slower using the trigonometric polynomials 
than using the basis functions. Thus, even for a basin 
with simplified geometry (such as a rectangle) basis 
functions are still preferred. 

1) Determination of the optimal mode 
number K opt 

The optimal mode number K opt strongly depends on 
the sampling strategy and the level of noise in the ob¬ 
servation. To illustrate this, the day 5 observational data 
are acquired from disposition B (Fig. 13b) with a red 
noise ( r] l = 0.5, R cor = 70 km) and the day 90 obser¬ 
vational data are acquired from disposition A (Fig. 13a) 
with a white noise (r] l = 0.5). Using VCCF, the optimal 
truncation K opt is determined as 25-30 (50 to 75) for 
day 5 (day 90) circulation. 

Taking the regularization parameter 8 [see (4.15) in 
Part I] as 0.1, the noisy data are reconstructed on day 
5 (Fig. 14a) and day 90 (Fig. 14c). Comparison between 
Figs. 14a and 12a (Figs. 14c and 12b) shows that the 
reconstruction scheme reproduces almost all the details 
of circulation pattern. The residual velocity vector fields 
for day 5 and day 90 are represented in Figs. 14b and 
14d, respectively. The reconstruction scheme underpre¬ 
dicts the basin wide circulation for day 5. 

2) Effect of noise on the reconstruction 
skill 

To test the sensitivity of reconstruction skill on the 
noise level and characteristics, six different types of 
noise are added to the exact observational data with 
various noise levels (0 < r\ l < 1.0) on day 90: (a) white 
noise, (b) color noise with R cor = 7.5 km, (c) color noise 
with R cor = 70 km, (d) color noise with R cor =100 km, 
(e) color noise with R cor = 200 km, and (f) color noise 
with R cor = 500 km. The velocity field is reconstructed 
using the schemes presented in Part I, and the rms errors 
of u component (Fig. 15) and v component (Fig. 16) 
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Fig. 14. (a) Velocity field and (b) residue velocity on day 5 reconstructed (8 = 0.1 and K opt = 30) from observed data collected from disposition B with a color noise r] l = 0.5, R COI = 
km; (c) velocity field and (c) residue velocity on day 90 reconstructed (8 = 0.1 and K opt = 50) from observed data collected from disposition A with a white noise with rj 1 = 0.5. 
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Fig. 15. Dependence of rms errors for u component (on day 90) 
on rj 1 for different noise characteristics: (a) white noise, (b) color 
noise with R COI = 7.5 km, (c) color noise with R COI = 70 km, (d) color 
noise with R cor =100 km, (e) color noise with R COI = 200 km, and 
(f) color noise with R COI = 500 km. 

are computed for different noise characteristics. The re¬ 
sults show that the reconstruction skill is quite high (x u 
< 0.3 and Xv < 0.2) and weakly sensitive to and to 
correlation radius of color noise if < 0.5 and R cor < 
500 km (Figs. 15 and 16). 

3) Normal velocity at the open boundary 

To verify the reconstruction quality at the open 
boundary, the reconstructed (denoted by +) and exact 
normal (solid curve) velocities at the open boundary are 
presented for various noise characteristics (noise level 
and correlation scale for color noise) and sampling strat¬ 
egies: disposition A (Fig. 17), disposition B (Fig. 18), 
and disposition D (Fig. 19). Among the six patterns (Fig. 
13), dispositions A, B, and C are domain wide sampling 
(some observations close to the open boundary); dis¬ 
position D is a partial domain sampling; and dispositions 
E and F are coastal sampling (all observations far away 
from the open boundary). 

( i) Effect of sampling 

For domainwide sampling (Figs. 17 and 18), the nor¬ 
mal velocity is reproduced quite well for evident noises 
( r] l = 0.5). However, for partial domain sampling (Fig. 
19), the normal velocity is reproduced only for the dis¬ 
position D with weak ( r\ x = 0.1) white noise, and not 
for the coastal sampling (dispositions E and F) where 
the observations are far from the open boundary. Such 
a failure is due to the fact that the developed noise 
reduction procedure (3.6)-(3.8) cannot totally eliminate 
the noise caused by the truncation error. 

Generally, the reconstruction skill does not directly 
depend on the number of observations, but does depend 
on the conditional number of the coefficient matrix A 



Fig. 16. Same as in Fig. 15, except for v component. 


of the linear algebraic equation [(4.5) in Part I]. This 
number is determined by the station disposition pattern. 
The farther the observations are to the open boundary, 
the larger the conditional number will be. It reaches 10 8 
for dispositions E and F. For a high conditional number, 
a weak noise in the observation will destroy the recon¬ 
struction. To improve this, the approach should be com¬ 
bined with a hydrodynamic model, which is the ongoing 
project. 

(ii) Effect of deterministic disturbance for 
domainwide sampling 

Without stochastic noise (Figs. 17 a,b; Fig. 18a), the 
difference between reconstructed and exact normal ve¬ 
locities shows the effect of a deterministic disturbance 
(i.e., the effect of the mode truncation) is caused by the 
truncation error. For observations sampled on day 90 
with disposition A (Fig. 17), the reconstruction skill is 
almost the same between a mode truncation of 50 (Fig. 
17a) and 70 (Fig. 17b), which indicates that the recon¬ 
struction skill is not sensitive to K opt within the range 
(50-75 for day 90) determined using VCCF. 

With a color noise (= 0.5, R cor = 70 km) in ob¬ 
servations sampled on day 5 with disposition B (Figs. 
18b,d) the reconstruction skill without regularization (8 
= 0) is lower with K opt = 21 (Fig. 18d) than with K opt 
= 30, which indicates that the reconstruction skill is 
sensitive to K opt when it varies out of the range (25-30 
for day 5) determined using VCCF. 

(iii) Effect of stochastic noise characteristics for 
domainwide sampling 

The reconstruction skill is not sensitive to noise char¬ 
acteristics if the mode truncation K opt varies within the 
range (25-30 for day 5, and 50-75 for day 90), deter¬ 
mined by VCCF (all the panels in Figs. 17 and 18 except 
Fig. 18d). The small difference between Figs. 17c,e,g, 
and Figs. 17d,f,h shows that the regularization (§ ^ 0) 
does not significantly improve the reconstruction skill 
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Fig. 17. Reconstructed (from disposition A, denoted by “ + ”) and exact (solid curve) normal velocities at the open boundary on day 90 
(a) without stochastic noise K opt = 50; (b) without stochastic noise K opt = 70; (c) with a white noise rj l = 0.5, and K opt = 56, 8 = 0; (d) 
with a white noise r] 1 = 0.5, and K opt = 56, 8 = 0.1; (e) with a color noise r\ x = 0.5, R COI = 100 km, and K opt = 50, 8 = 0; (f) with a color 
noise r] 1 = 0.5, R COI = 100 km, and K opt = 50, 8 = 0.1; (g) with a color noise r} l = 0.5, R cor = 200 km, and K opt = 50, 8 = 0; (h) with a 
color noise, rj 1 = 0.5, R cor = 200 km, and K opt = 50, 8 = 0.1. 


without regularization (5 = 0) for domain wide sampling 
with high density. Notice that the mode truncation K opt 
in all cases illustrated in Fig. 17 is within the range (50- 
75 for day 90) determined using VCCF. 

4. Case 3: Reconstruction from Lagrangian drifter 
data in the Southern Ocean 

The large-scale circulation is reconstructed from the 
observations of more than 200 surface drifting buoys 
deployed during FGGE between 1 December 1978 and 
1 December 1979 to illustrate the capability of the 
scheme on processing the Lagrangian-type data. 

a. FGGE data 

The original buoy data consist of nonuniform time 
series of position, sea level pressure, and sea surface 


temperature. Garrett (1980) discussed the accuracy of 
the measurements and buoy data acquisition system. The 
long-term trajectories of these buoys in the Southern 
Ocean and surrounding waters (south of 25°S) are ex¬ 
amined (Fig. 20). 

The surface velocity from a drifter buoy V drift can be 
represented by 

^drif't ^geostr ^Ekm ^Stok ^drag’ 1 ) 

where V geostr , V Ekm , V Stok , and V drag are the geostrophic, 
Ekman, Stokes, and drag velocities, respectively. Dan- 
iault and Menard (1985) estimated that the mean surface 
velocity (V drift ) over the Southern Ocean is about 24 cm 
s -1 ; and the summation of the Ekman, Stokes, and drag 
velocities (i.e., the ageostrophic velocity) is around 7 
cm s -1 . Therefore, the ageostrophic component is sig¬ 
nificantly large (around 30%) in the drifter buoy data. 
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Fig. 18. Reconstructed (from disposition B, denoted by “ + ”) and exact (solid curve) normal velocities 
at the open boundary on day 5 (a) without stochastic noise K opl = 30; (b) with a color noise r] 1 = 0.5, R COI 
= 70 km, and K opt = 30, 8 = 0; (c) with a color noise r] l = 0.5, R COT = 70 km, and K opt = 30, 8 = 0.1; (d) 
with a color noise r] l = 0.5, R COT = 70 km, and K opt = 21, 8 = 0. 


Thus, both scalar potentials T' and <F are used to re¬ 
construct the data. 

b. Data preparation 

The trajectories of the drifter buoys were smoothed 
using 2D splines to exclude any spatial peculiarities 
along the trajectories, which cannot be resolved on a 2° 
X 2° grid. Figure 21 shows the Lagrangian velocities 
along the smoothed trajectories for May 1979. After the 
smoothing, the ratio between the filtered ( E mu ) and total 
(E tot ) kinetic energies is given by 

— - 0.2 - 0.25. 

E tot 


The domain B (Fig. 22b) is enclosed by two rigid 
boundaries: T and T 3 (coastline of South America south 
of 48°S) and an open boundary T' (48°S latitude over 
water). Comparison between the two reconstructed ve¬ 
locity fields leads to the verification of the open bound¬ 
ary condition [(2.10) of Part I] and the boundary ex¬ 
tension method (section 4e of Part I). 

d. Spectral decomposition 

Both domains are multiply connected. Here, domain 
A has two islands (enclosed by T and TJ, and domain 
B has one island (enclosed by T), respectively. Following 
section 2d(3) of Part I, the toroidal and poloidal spectral 
expansions for domains A and B are represented by 


c. Domains A and B 


\LTA — n \|/(1) —L p \Lr(2) + V n \Lta 

1 harm 1 ^2 * harm 1 /Li u k^k'> 


We reconstruct the circulation separately in two areas: 
domains A and B. The domain A (Fig. 22a) is enclosed 
by three rigid boundaries: T (70°S); T 2 (25°S and con¬ 
nected coastlines of South America, South Africa, Mad- 
cascar, and Australia south of 25°S); and T x (repre¬ 
senting New Zealand). Notice that the boundary con¬ 
ditions for domain A are often used in large-scale ocean 
circulation modeling for the Southern Ocean. 


v* = + 2 a**?, 


(4.2) 


d<P A 

dz 


V b <P A , 

/ -t m mi 


m —0 


dz 


y b 

/ ^ m m ’ 


m=0 


(4.3) 


where (T^^, T^f) and (<F A , O^) are the eigenfunctions of 
the Laplacian operator calculated for the two domains 
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Fig. 19. Reconstructed (from disposition D, denoted by “ + ”) and exact (solid curve) normal velocities at 
the open boundary white a noise (r] 1 = 0.1) on (a) day 5 with K opl = 28, and (b) day 90 with ^ opt = 75. 


with homogeneous boundary conditions [(2.24), (2.27), 
and (2.28) in Part I]. The additional three harmonic 


= 0, 

mjr, = 

i, 

¥fiLlr 2 = 0, 

functions ''LA,,, and satisfy the Laplacian 

*£Llr 

= 0, 

= 

i, 

^Ur, = 0, 

equation with the boundary conditions written by 

^hamJr 

= o, 

^h2Jr 3 ur' = 

i. 

(4.4) 



Fig. 20. Original FGGE buoy trajectories in the Southern Ocean. 


e. Determination of coefficients 

The spectral coefficients C l9 C 2 , and C 3 can be de¬ 
termined from data or using special integral relation¬ 
ships (Kamenkovick 1961; McWilliams 1977; Flierl 
1977; Platzman 1979) around any island: 



a = r, r„ r 2 , or r 3 u r, (4.5) 

where F T is the circulation along any island boundary. 
It is difficult to use (4.5) if F T is unknown. 

The constants C l5 C 2 , and C 3 are directly calculated 
from the drifter data together with the spectral coeffi¬ 
cients in (4.2). Notice that (4.5) can be easily included 
into the reconstruction process using Lagrangian mul¬ 
tipliers. 

f Regression equations 

For domain A, the linear regression is used in the 
reconstruction. For domain B, the nonlinear regression 
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Fig. 21. Lagrangian velocities along the smoothed buoy 
trajectories in May 1979. 


model [(4.2)-(4.5) of Part I] is used since the basis 
functions {T^} depend on the constant C 3 . 

g. Reconstructed circulation in domain A 

Domain A is enclosed by rigid boundaries. The basis 
functions of both potentials are obtained using the rigid 
boundary conditions [(2.7) of Part I]. For domain A, 
1596 observations are analyzed from 213 buoys, and 
the VCCF is used to determine the model parameters 
(K opt , M opt ) by (30, 4). The reconstructed circulation 
(Fig. 23) shows the main characteristics of the Southern 
Ocean circulation, such as the Antarctic Circumpolar 
Current, the Brazil-Malvinas Confluence, the Agulhas 
Current, the West Australian Current, as well as the 
South Pacific gyre, the Indian Ocean gyre, and the South 
Atlantic gyre. The Antarctic Circumpolar Current me¬ 
anders and has several branches extending into the Pa¬ 
cific Ocean. The mean current speed over the Southern 
Ocean is about 20 cm s _1 . Since the data are only on 
the sea surface, it is not estimated the volume transport 
in the Southern Ocean, such as through the Drake Pas¬ 
sage. 

h. Reconstructed circulation in domain B 

Domain B is a subarea of domain A, and has an open 
boundary along 48°S. The iteration process [(4.2)-(4.4) 
of Part I] are used to get the spectral coefficients. The 
boundary extension method is used to obtain the first 
guess of the spectral coefficients 8h m for the iteration. 
The calculated velocity along the open boundary (48°S) 
is used to determine the coefficient k for this domain. 




Fig. 22. (a) Domain A with two rigid boundaries T, T 2 , and an 
island C; and (b) domain B with two rigid boundaries T, r 3 , and 
one open boundary f\ 


For the reconstruction of domain B, 638 observations 
from 112 buoys are analyzed; VCCF is used to deter¬ 
mine the optimal mode truncation (K opt , M opt , S opt ) as 
(30, 4, 7); and a seventh-order polynomial is used to 
smooth k along the open boundary. 

Thus, two velocity fields for domain B are obtained: 

(a) a subset of the domain A velocity field (Fig. 24a), and 

(b) a reconstructed field for domain B, only with an open 
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boundary at 48°S (Fig. 24b). The two flow fields are quite 
similar. The residue (Fig. 25) between the two fields (sub¬ 
set of domain A minus domain B) is quite weak, and the 
relative root-mean-square difference between the two 
fields is calculated using (5.10) of Part I: 

- 0.09 and - 0.13, (4.7) 

which indicates that the scheme leads to 9% to 13% of 
errors in processing the drifter buoy data. 

The reconstructed normal and tangential velocities 
along 48°S (open boundary for domain B) obtained from 
domain A’s reconstruction (denoted by “1”) using the 
rigid boundary conditions (for a larger domain) and do¬ 
main B’s reconstruction (denoted by “2”) using the it¬ 
eration procedure [(4.2)-(4.4) of Part I] are examined 
(Fig. 26). The difference between the two fields for both 
components is quite small, except for the tangential 
component in the South Pacific sector (150°E-100°W). 
The maximum difference is found at 165°W: the tan¬ 
gential velocity is 25.5 cm s _1 from domain As esti¬ 
mation, and about 13 cm s _1 from domain B’s esti¬ 
mation. 


5. Conclusions 

First, the results from three experiments are presented 
to illustrate the capability of the theoretical approach 
developed in Part I in flow reconstruction from noisy 
and sparse data: (a) reconstruction of three-dimensional 
velocity at 20-m depth in the eastern Black Sea, (b) 
reconstruction of circulation pattern in a wind driven 
shallow water basin flow with open boundary, and (c) 
reconstruction of large-scale circulation in the Southern 
Ocean from the FGGE drifter data. 




Fig. 24. Reconstructed circulation in domain B (a) obtained from 
Fig. 23 and (b) calculated using the open boundary condition (4.10) 
of Part I. 


Second, the difficult problem in any reconstruction is 
how to select the optimal field from noisy data. In ocean¬ 
ography, the low-order noise statistics are usually un¬ 
known. This difficulty is overcome using a special reg¬ 
ularization algorithm. The three cases presented dem¬ 
onstrate the capability of the algorithm in reconstruction 
with the noise-to-signal ratio and the condition number 
up to 1 and 10 8 to 10 10 , respectively. 

Third, the results from the three cases demonstrate 
that the sampling strategy is the key issue for the re¬ 
construction skill. In the eastern Black Sea experiment, 
the reconstruction skill deteriorates with the reduction 
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Fig. 25. Residual velocity vectors between the two calculations. 


of sampling density from high sampling density (19 
observations) to medium sampling density (14 obser¬ 
vations), and from medium sampling density to low 
sampling density (9 observations). For the low sampling 
density, we can only estimate six toroidal modes and 
one poloidal mode. Any combination of these modes 
results to a poor reconstruction skill. 

Fourth, the empirical functional / emp is not the best 
cost function for estimating the optimal mode truncation 
of the toroidal and poloidal potentials or other free ex¬ 
ternal model parameter from limited data. It is recom¬ 
mended to use the Vapnik-Chervonkis approach (Vap- 
nik 1983) instead of / emp since the VCCF is proven in 
the three cases to be robust to the variations of the 
sampling size. 

Fifth, reconstruction skill depends on free model pa¬ 
rameters (K opt , M opt , S opt ). The three cases show that the 
reconstruction of the normal velocity at the open bound¬ 
ary is more accurate than the tangential velocity. How¬ 
ever, this result might not be universal since the recon¬ 
struction skill of normal and tangential velocities de¬ 
pends on the flow structure along the open boundary. 

Sixth, the reconstruction process is reduced to a linear 
regression model for known and unknown simplified 
boundary conditions if the calculation is over a simply 
connected domain. For a multiply connected domain 
with unknown or known boundary conditions, the flow 
reconstruction is only possible using a nonlinear re¬ 
gression model with independent parameterization for 
boundary conditions. 
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